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Abstract 

In an attempt to explain the latest experimental result about the branching 
ratio of to — ► 7r°7r°7 decay we reexamine the decay mechanism of this decay 
in a phenomenological framework in which the contributions of VMD, chiral 
loops, a- meson intermediate state amplitudes and the effects of to — p mixing 
are considered. We conclude that in order to obtain the experimental value 
of the branching ratio B(u — > 7r 7r°7) cr-meson amplitude which makes a 
substantial contribution should be included into the reaction mechanism and 
the effects of u — p mixing should be taken into account. We also estimate the 
coupling constant g WCJ1 as g wal = 0.11 which is much smaller than the values 
suggested by light cone QCD sum rules calculations. 
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The recent experimental study of p — > 7r 7r°7 and uj — > 7r°7r°7 decays by SND Collabo- 
ration obtained the value B(u — > 7r°7r°7) = (6.6^o;§ ± 0.6) x 10~ 5 for the branching ratio 
of the uj — > 7r°7r°7 decay [1]. Their result is in good agreement with GAMS Collaboration 
measurement of B(uo — > 7r°7r°7) = (7.2 ± 2.5) x 10~ 5 [2], but it has a higher accuracy. 

On the theoretical side, uj — > 7r°7r°7 decay was first studied by Singer [3] who postulated 
that this transition proceeds through the uj — > (p)7r — > 7r 7r°7 mechanism involving p-meson 
intermediate state. The contribution of intermediate vector mesons (VMD) to the vector 
meson decays into two pseudoscalars and a single photon V — > PP'j was also considered by 
Bramon et al. [4] using standard Lagrangians obeying the SU(3)-symmetry, and in particular 
for the branching ratio of the decay uj — > 7r°7r°7 they obtained the result B(uo — > 7r°7r°7) = 
2.8 x 1CT 5 . The V — > PP lr y decays have also been considered within the framework of chiral 
effective Lagrangians using chiral perturbation theory. Bramon et al. [5] studied various such 
decays using this approach and they noted that if chiral perturbation theory Lagrangians are 
used there is no tree-level contribution to the amplitudes for the decay processes V — > PP'^f, 
and moreover the one-loop contributions are finite and to this order no counterterms are 
required. They considered both tttt and KK intermediate loops. In the good isospin limit 
7r-loop contributions to uj — > 7r°7r°7 amplitude vanish, and the contribution of K-loops is 
two orders of magnitude smaller than the contribution of VMD amplitude. Therefore, the 
VMD amplitude essentially accounts for the decay rate of the uj — > 7r°7r°7 decay. Guetta 
and Singer [6] recently updated the theoretical value for the branching ratio B(uo — > 7r°7r°7) 
of the decay uj — > 7r°7r°7 as B{uj — > 7r°7r°7) = (4.1 ± 1.1) x 1CT 5 . In their calculation they 
noted that when the Born amplitude for VMD mechanism is used the decay rate uj — > 7r°7r°7 
is proportional to the coupling constants and <?L. 7 , and they assumed that the decay 
uj — > 3n proceeds with the same mechanism as uj — > 7r°7r°7, that is as uj — > (p)n — > nnn 
[7]. They use the experimental inputs for the decay rates T(uj — * 3ir), T(p° — > 7r°7) and 
r(p — > 7T7r), and furthermore they employ a momentum dependent width for p- meson. If 
a constant p-meson width is used, then the value B(uj — > 7r°7r°7) = (3.6 ± 0.9) x 10~ 5 is 
obtained for the branching ratio of this decay. Therefore there appears to be a serious 
discrepancy between the theoretical result and the experimental value for the branching 
ratio of the uj — > 7r°7r°7 decay. 

Guetta and Singer [6] noted that in the theoretical framework based on chiral pertur- 
bation theory and vector meson dominance one feature has been neglected. This is the 
possibility of uj — p mixing one consequence of which is the isospin violating uj — > 7r + 7r~ 
decay with the branching ratio B(uj — > 7r + 7r~) = (2.21 ± 0.30)% [8]. The phenomenon of 
uj — p mixing has been observed in the electromagnetic form factor of the pion improving the 
standard VMD model result involving p-meson intermediate state. A recent review of uj — p 
mixing and vector meson dominance is given by O'Connell et al. [9]. Guetta and Singer [6] 
calculated the effect of uj — p mixing using the Born amplitude for VMD mechanism and 
showed that it increases the uj — > 7r°7r°7 width only by 5%. They then combined all the 
improvements on the simple Born term of VMD mechanism, that is uj — p mixing, momen- 
tum dependence of p-meson width and the inclusion of the chiral loop amplitude as given 
by Bramon et al. [5] , and using the resulting amplitude which includes all these effects they 
obtain the theoretical result B(uj — > 7r°7r°7) = (4.6 ± 1.1) x 10 -5 for the branching ratio of 
the uj — > 7r°7T 7 decay. Palomar et al. [10] also analyzed the radiative V — > PP'j decays 
using vector meson dominance, chiral loops obtained using unitarized chiral perturbation 
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theory, and uj — p mixing. They obtained the result B{uj — > 7r°7r°7) = (4.7 ± 0.9) x 10~ 5 for 
this branching ratio. These theoretical results are still seriously less than the latest experi- 
mental result B(uj — > 7r°7r°7) = (6.6tJ;g±0.6) x 10~ 5 . Therefore, the possibility of additional 
contributions to the mechanism of uj — > 7r°7r°7 decay should be investigated. 

One such additional contribution to the uj — > 7r°7r°7 decay may be provided by the 
amplitude involving scalar-isoscalar er-meson as an intermediate state. The existence of the 
controvertial cx-meson now seems to be established by the Fermilab (E791) Collaboration in 
their observation of the D + — > <jtt + — > tttttt decay channel in which cx-meson is seen as a clear 
dominant peak with M a = (438 ±31) MeV, and T a = (338 ±48) MeV [11], where statistical 
and systematic errors are added in quadrature [12]. Two of present authors in a previous 
work [13], calculated the decay rate for the decay uj — > 7r°7r°7 by considering p-pole vector 
meson dominance amplitude as well as cx-pole amplitude in a phenomenological approach. 
By employing then the available experimental value for this branching ratio B(uj — > 7r°7r°7) = 
(7.2 ± 2.5) x 1CT 5 [2] which is somewhat less accurate than the present new value [1], they 
obtained for the coupling constant g W(n the values g W(ri = 0.13 and g W(77 = —0.27. They 
observed that cx-meson intermediate state amplitude makes an important contribution by 
itself and by its interference with the VMD amplitude. The same authors later estimated 
the coupling constant g ual by studying cjcx7-vertex in the light cone QCD sum rules method 
[14], and the value | g wery \— (0.72 ± 0.08) was deduced for this coupling constant. Aliev et 
al. [15] also used light cone QCD sum rules techniques to calculate the coupling constant 
g pal , and they obtained the value g pal = (2.2 ±0.2) from which by using SU(3)-symmetry it 
follows that the coupling constant g U(T1 should have the value g wai = 0.73. Thus, there seems 
to be a serious discrepancy between the values obtained for the coupling constant g wr7 ^ using 
light cone QCD sum rules method and the phenomenological analysis of the uj — > 7r°7r°7 
decay. 

Therefore, in the present work, we reconsider the uj — > 7r°7r°7 decay in a phenomenological 
framework in order to the assess the role of cx-meson in the mechanism of uj — > 7r°7r°7 decay 
and to recalculate the coupling constant g^ cl utilizing the latest experimental value of the 
branching ratio uj — > 7r°7r°7. For this purpose, we calculate the decay rate for the decay 
uj — > 7r°7r°7 by considering p-pole vector meson dominance amplitude, chiral loop amplitude, 
cr-pole amplitude and we also include the effects of uj — p mixing which was not taken into 
account in the previous analysis [13]. 

In order to calculate the effects of the uj — p mixing in the uj — > 7r°7r°7 decay we need 
an amplitude characterizing the contribution of cx-meson to the p° — > 7r°7r°7 decay. In 
a previous work, two of the present authors [16] calculated the branching ratio B(p° — > 
7r + 7r~7) in a phenomenological framework using pion bremsstrahlung amplitude and cx-meson 
pole amplitude, and they determined the coupling constant g pul by using the experimental 
value of the branching ratio B(p° — > 7r + 7r~7). In a following work, these authors [17] 
calculated the branching ratio B(p° — > 7r°7r°7) of the p° — > 7r°7r°7 decay using the value 
of the coupling constant g pai they thus obtained again in a phenomenological approach in 
which the contribution of cx-meson, w-meson intermediate states and of the pion-loops are 
considered. However, the branching ratio B(p° — > 7r°7r°7) obtained this way was much 
larger than the experimental value. This unrealistic result was due to the constant p° — > CX7 
amplitude employed and consequently the large coupling constant g pal that was deduced 
using the experimental branching ratio of the p° — > 7r + 7r~7 decay. Therefore, it may be 
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concluded that it is not realistic to include a-meson in the mechanism of radiative p°-meson 
decays as an intermediate pole state. On the other hand, the V — > PP lr j decays were also 
studied by Marco et al. [18] in the framework of chiral perturbation theory. They used the 
techniques of chiral unitary theory [19], and by a unitary resummation of the pion loops 
through the Bethe-Salpeter equation they obtained the decay rates for various decays. They 
furthermore noted that their result for p° — > 7r°7r°7 decay could be interpreted as resulting 
from the p° — > ((7)7 — > 7r°7r°7 mechanism. Thus, it seems that a natural way to include the 
effects of cr-meson in the mechanism of radiative p°-meson decays is to assume that cr-meson 
couples to p°-meson through the pion-loop. 

Our phenomenological approach is based on the Feynman diagrams shown in Fig. 1 for 
uj — > 7r°7r°7 decay and in Fig. 2 for p° — > 7r°7r°7 decay. The direct terms shown in the 
diagrams in Fig. 1 b and in Fig. 2 b, c are required to establish the gauge invariance. The 
interaction term for two vector mesons and one pseudoscalar meson is given by the Wess- 
Zumino anomaly term of the chiral Lagrangian [20], we therefore describe the ujp-n- vertex 
by the effective Lagrangian 

£1% = 9 Ufm ^ a %^d a p p ■ 7? , (1) 

which also conventionally defines the coupling constant g Ufm . This coupling constant was 
determined by Achasov et al. [21] through an experimental analysis as g Ufm = (14.4 ± 
0.2) GeV~ 1 , who assumed that uj — > 37r decay proceeds with the intermediate pir state as 
uj — > (p)7r — > nnn and they used the experimental value of the uj — > 3n width to deduce the 
coupling constant g paJ7r . Similarly we describe the Vn^- vertices where V = p, uj with the 
effective Lagrangian 

C^l r = g Vin €^ a %V v d a Apn . (2) 

We then use the experimental partial widths of the radiative V — > 717 decays [8] to deduce 
the coupling constants g umi and g p7ry - This way for the coupling constants g W7r7 and g p7T1 we 
obtaine the values g^ = (0.706 ± 0.021) GeV' 1 and g p7TJ = (0.274 ± 0.035) GeV' 1 . The 
p7T7r-vertex is described by the effective Lagrangian 

L% = W,-(^7fX7f) , (3) 

and the experimental decay width of the decay p — > nn [8] then yields the value g pw7T = 
(6.03 ± 0.02) for the coupling constant g fm7T . We describe the <77T7r-vertex by the effective 
Lagrangian 

CJI = \g^M^ ■ ira . (4) 

The experimental values for M a and T a M a = (483 ±31) MeV and T a = (338 ± 48) MeV, 
where statistical and systematic errors are added in quadrature [11,12] then results in the 
strong coupling constant g an7T = (5.3 ± 0.55). The effective Lagrangians C e JJ % and £ e p H 
are obtained from an extension of the a model where the isovector p is included through 
a Yang-Mills local gauge theory based on isospin with the vector meson mass generated 
through the Higgs mechanism [22] . In order to describe the 7r 4 - vertex we again consider the 
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cr-model with spontaneous symmetry breaking [23] and we describe the 7r 4 -vertex by the 
effective Lagrangian 



C eff 



t(") 2 



where the coupling constant A is given as A = — 9n ^ N MrT M ^ 7r and we use 97rNN 



Air 



(5) 



14. We note 



that this effective interaction results in only isospin 1=0 amplitudes. The small 1=2 am- 
plitudes were also neglected in previous calculations within the framework of chiral unitary 
theory [18]. The uoa^- vertex is described by the effective Lagrangian 

C e £ r = ^-9 u<n ff x u^d a Ap-dpA a )a , (6) 

which also defines the coupling constant <3W 7 . 

In our calculation of the invariant amplitude, we make the replacement q 2 — M 2 — > 
q 2 — M 2 + iMT in p-meson and <7-meson propagators. We use for cr-meson the momentum 
dependent width that follows from Eq. 4 

, M 2 



° q 2 \M 2 - AM 2 



AM 2 



(7) 



and for p-meson we use the following momentum dependent width as conventionally adopted 
[9] 
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Loop integrals similar to the ones appearing in Figs. 1 and 2 were evaluated by Lucio 
and Pestiau [24] using dimensional regularization. We use their results and, for example, 
we express the contribution of the pion-loop amplitude corresponding to p° — > (7r + 7r~)7 — > 
7r°7r°7 reaction in Fig. 2 b as 

6 QpTVTV^ 



2tt 2 M 2 
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where a = M 2 /M 2 , b = (p — k) 2 /M 2 , p(u) and k(e) being the momentum (polarization 
vector) of p-meson and photon, respectively. A similar amplitude corresponding to p° — > 
{it + -K~)^a — > 77r°7r° reaction can also be written. The function I(a,b) is given as 
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Although the contribution of chiral kaon-loop diagram shown in Fig. 1 b to the decay rate 
of uo — > 7r°7r°7 decay is small, we also include the corresponding amplitude of this diagram 
in our calculation for completeness. However, since we lack any experimental information 
to describe the uoK + K~ -vertex and K + K~-n°ir amplitude, for the contribution of this 
diagram we use the amplitude given by Bramon et al. [5] derived using chiral perturbation 
theory. This may not be entirely consistent with the philosophy of our phenomenological 
approach, but since its contribution is shown to be small [5] we do not think that this 
way of including kaon-loop diagram into our calculation constitutes a serious inconsistency. 
Moreover, in Fig. 2 b in addition to pion-loop intermediate state there is also a contribution 
to p° — > 7r°7r°7 decay coming from KK intermediate state. However, as shown by Bramon et 
al. [5], these kaon-loop intermediate states give a contribution to p° — > 7r°7r°7 decay which is 
10 3 times smaller than the contribution coming from the charged-pion loops. Therefore, in 
our calculation we do not take the kaon-loop amplitude in p° — > 7r°7r°7 decay into account. 

We describe the uo — p mixing by an effective Lagrangian of the form 

Cf w = Yl pLO uo^ , (12) 

where uo^ and p p denote pure isospin field combinations. The corresponding physical states 
can therefore be written as [9] 

I p > = \p,I=l>+e\uj,I = 0> 

\ uo > = \ uo, I = > -e \ p, I = 1 > , (13) 



where 



m*-w • /.u„r„ iMJ', 



(14) 



O'Connell et al. [9] determined U puJ from fits to e + e~ — > 7r + 7r~ data as Yi pu) = (—3800 ± 
370) MeV 2 . Then, using the experimental values for My and T v , the mixing parameter e 
is obtained as e = (—0.006 + «0.036). Another effect of uo — p mixing besides the mixing of 
the states is that it modifies the p-propagator in diagrams in Fig. 1 a as 
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D p (s) D p (s) \ g p7rl D p (s) 



(15) 



where D p (s) = s — M p + iM p T p (s). The amplitude of the decay lj — > n ir 7 can then be 
written as M. = M.® + eA4' where A4q includes the contributions coming from the diagrams 
shown in Fig. 1 for lj — > 7r°7r°7 and M' represents the contributions of the diagrams in Fig. 
2 for p° -> 7r°7r°7. 

We calculate the invariant amplitude A / i(E 7 , Ei) this way for the decay uo — > 7r°7r°7 from 
the corresponding Feynman diagrams shown in Fig. 1 and 2 for the decays uo — > 7r°7r°7 and 
p° — > 7r°7r°7, respectively. The differential decay probability for an unpolarized u;-meson at 
rest is then given as 

M | 2 , (16) 



dE 1 dE 1 (2vr) 3 8M W 
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where E 7 and Ei are the photon and pion energies respectively. We perform an average over 
the spin states of the vector meson and a sum over the polarization states of the photon. 
The decay width is then obtained by integration 

1 \ / 7 ,max. rEi ,max. (J I 

where the factor (|) is included because of the 7r°7r° in the final state. The minimum 
photon energy is E 7>m j n . = and the maximum photon energy is given as E 1 ^ max- = (M% — 
4M^)/2M Ct) =341 MeV. The maximum and minimum values for pion energy Ei are given by 

1 

n2 



[-2E*M U + 3E,M^ - Mi 



2(2/-;- \L-\n\ 

±E^(-2E,M LU + M^)(-2E 1 M U) + Ml - AMI) ] . 

The theoretical decay rate for uj — > 7r°7r°7 decay that we calculate results in a quadric 
equation for the coupling constant g U(n , and using the experimental value for this decay rate 
[1], we obtain the values g wa ^ = (0.11 ± 0.01) and g wo - 7 = (—0.21 ± 0.02) for the coupling 
constant g u<n . These values are smaller than the values g u<n = 0.13 and g W(J1 = —0.27 that 
were obtained in the previous phenomenological analysis which did not include the effects 
of uj — p mixing [12]. We may therefore conclude that uj — p mixing does make a reasonably 
substantial contribution to the uj — > 7r°7r°7 decay amplitude when a-meson intermediate 
state is taken into account which then results in a reduced value for the coupling constant 

guja-y- 

The resulting 7r°7r°-invariant mass distribution for the uj — > 7r°7r°7 decay if we use the 
coupling constant g uaj = 0.11 is plotted in Fig. 3 where we also indicate the contributions 
coming from the different amplitudes. The interference term between the different ampli- 
tudes is positive over the whole region. The contribution of the VMD amplitude calculated 
from the diagrams in Fig. 1 a does not change appreciably if we take into account the 
effect of uj — p mixing by including the contribution coming from the diagrams in Fig. 2 a 
The situation changes somewhat if we consider VMD and chiral loop amplitudes and uj — p 
mixing as well. However, the significant modification is obtained when we include VMD, 
chiral loop, and cr-meson intermediate state amplitudes with uj — p mixing. 

Although the shapes of the various curves are quite similar, as it can also be seen the cor- 
responding branching ratios are considerably different. If we use only the VMD amplitudes 
we obtain for the branching ratio B(uj — > 7r°7r°7) the values B(uj — > 7r°7r°7) = 3.96 x 10~ 5 and 
B{uj — > 7r°7r°7) = 4.22 x 10~ 5 without and with the effect of the uj — p mixing, respectively. 
These results are in agreement with the calculation of Guetta and Singer [6]. When we 
consider VMD amplitudes and chiral loop amplitudes, the resulting values for the branching 
ratio are B{uj -> n°ir°~f) = 3.98 x 10~ 5 and B{uj -> 7r°7r°7) = 4.67 x 10~ 5 without and with 
the effect of the uj — p mixing included, respectively. Again these calculated values are consis- 
tent with the previous results, in particular, with the results of Bramon et al. [5] and Palomar 
et al. [10]. However, the main difference with previous results arises when we consider the 
contribution of cr-meson intermediate state. Indeed, the branching ratio B(uj — > 7r°7r°7) that 
we calculate using the full amplitude including the contributions of VMD, chiral loop, and 
cr-meson intermediate state diagrams which is B{uj — > 7r°7r°7) = 7.29 x 10~ 5 is reduced to 
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the value B(oj — > 7r°7r°7) = 6.6 x 1CT 5 when the effect of the uj — p mixing is included. We 
thus observe that in this case the effect of uj — p mixing on the amplitude of the uj — > 7r°7r°7 
decay is reasonably pronounced and moreover the contribution of the cr-meson intermediate 
state is quite substantial. We may therefore conclude that cr-meson intermediate state and 
uj — p mixing should be included in the analysis of uj — > 7r°7r°7 decay in order to explain the 
latest experimental result [1]. 

In Fig. 4, we plot the 7r 7r°-invariant mass distribution we obtain if we use the coupling 
constant g ueri = —0.21. In this case also the cr-meson intermediate state amplitude makes 
a very significant contribution, but the interference term resulting from the cr-meson and 
VMD and chiral loop amplitudes is negative over some regions of the spectrum. Moreover, 
the overall shape is quite different from the previous case. 

In our work, we consider a new approach to the mechanism of p° — > 7r°7r°7 decay. We 
assume that the cr-meson couples to p°-meson through the pion-loop and thus we neglect a 
direct po-7-vertex. The assumed mechanism of p° — > 7r°7r°7 decay is shown by the Feynman 
diagrams in Fig. 2. We can therefore calculate the branching ratio of p° — > 7r°7r°7 decay 
based on these Feynman diagrams. This way we obtain for the branching ratio of p° — > 7r°7r°7 
decay the value B(p° — > 7r°7r°7) = (3.81 ± 0.63) x 10~ 5 which is in a good agreement the 
latest experimental result B(p° — > 7r°7r°7) = (4.1^^0.3) x 10~ 5 [1]. Scalar cr-meson effects 
in radiative p°-meson decays are studied in a recent work by two of the present authors in 
detail [25]. In that work, the standard 7r + 7r~ — > 7r°7r° amplitude of chiral perturbation 
theory is used in the loop diagrams, and the value B(p° — > 7r°7r°7) = (4.95 ± 0.82) x 10~ 5 
is obtained for the branching ratio. If we use the same amplitude in the present work we 
then obtain the value B(uj — > 7r°7r°7) = 6.12 x 10~ 5 using the full amplitude including the 
contributions of VMD, chiral loop, and cr-meson intermediate state diagram as well as the 
effects of uj — p mixing, which is not very different from B(lu — > 7r°7r°7) = 6.6 x 10~ 5 that is 
obtained employing the effective Lagrangian given in Eq. 5 to describe the 7r 4 -vertex. 

An essential assumption of our work, therefore, is that there is no SU(3) vector meson- 
sigma-gamma vertex. Thus, the u;cr7-vertex cannot be related to the po"7-vertex. The ua / ~f- 
vertex that we use may be considered as representing the effective final state interactions 
in the 7T7r-channel. The small value of the coupling constant g W(J1 that we obtain leads 
to a change in the Born amplitude of the uj — > 7r°7r°7 decay which is of the same of the 
magnitude, as it is typical of final state interactions. As a matter of fact, Levy and Singer 
[26] in a study of the uj — > 7r°7r°7 decay using dispersion-theoretical approach shown that 
final state interactions resulting in a decay rate of the same order of magnitude as the one 
calculated from the Born term can be parameterized with the effective-pole approximation. 

As a final comment, we like to mention that our analysis suggests that the coupling 
constant g^ ai has actually a much smaller value than obtained by light cone QCD sum rules 
calculations. 
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FIG. 1. Feynman diagrams for the decay uj — ► 7r 7r°7 
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FIG. 2. Feynman diagrams for the decay p — > 7r 7r°7. 
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FIG. 3. The tt it invariant mass spectrum of the decay uj — > 7r 7r 7 for g wo - 7 = 0.11. The 
separate contributions resulting from the amplitudes of VMD; VMD and uj — p mixing; VMD, 
chiral loop, uj — p mixing; VMD, chiral loop, a-meson intermediate state, uj — p mixing are shown. 
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FIG. 4. The it it invariant mass spectrum of the decay lo — ► it tt 7 for g^a-y = —0.21. The 
separate contributions resulting from the amplitudes of VMD; VMD and uj — p mixing; VMD, 
chiral loop, uj — p mixing; VMD, chiral loop, a-meson intermediate state, uj — p mixing are shown. 
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